Doped Singlet-Pair Crystal in the Hubbard model on the checkerboard lattice 
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In the limit of large nearest-neighbor and on-site Coulomb repulsions, the Hubbard model on 
the planar pyrochlore lattice maps, near quarter-filling, onto a doped quantum fully packed loop 
model. The phase diagram exhibits at quarter filling a novel quantum state of matter, the Resonating 
Singlet-Pair Crystal, an insulating phase breaking lattice symmetry. Properties of a few doped holes 
or electrons are investigated. In contrast to the doped quantum antiferromagnet, phase separation is 
restricted to very small hopping leaving an extended "window" for superconducting pairing. However 
the later is more fragile for large hopping than in the case of the antiferromagnet. 

PACS numbers: 75.10.-b, 75.10.Jm, 75.40.Mg, 74.20.Mn, 71.10.Fd 
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Correlated fermions and bosons on frustrated lattices 
exhibit fascinating properties. For example, on the tri- 
angular lattice, hard-core bosons with nearest neighbor 
(NN) repulsion can realize a supersolid phase exhibiting 
both charge order and superfluidity [J. In the light of 
recent STM experiments suggestingthe coexistence of su- 
perconductivity and charge order it is interesting to 
ask whether a fermionic analog of a supersolid could be 
realized by doping a Mott insulator whose ground state 
spontaneously breaks lattice symmetries ? 

In this Letter we study correlated fermions on the 
checkerboard lattice, a two-dimensional array of corner- 
sharing tetrahedra. This is a two-dimensional analog of 
the pyrochlore lattice found in numerous spinel and py- 
rochlore materials [Sj. For a filling of exactly half an 
electron per site, extremely strong on-site and NN repul- 
sion select a macroscopically degenerate manifold of low- 
energy configurations fulfilling the so-called "ice rules" 
constraint, i.e. having exactly two particles per tetrahe- 
dron. These can be mapped onto the configuration space 
of a 6- vertex model Q . By associating each particle with 
a dimer joining the centers of the two corner-sharing 
tetrahedra, an alternative fully packed loop representa- 
tion can be obtained. Although the constrained quantum 
dynamics of bosons [3| and spinless fermions 01 differ, the 
phase diagrams of these models contain a rich variety of 
crystalline phases breaking lattice translations and/or ro- 
tations. In addition, fractionalization of a doped charge 
e into two e/2 components can also appear under some 
conditions [1, 0, Hi- 

Here, we consider the realistic case of spinful fermions 
(electrons) whose spin degrees of freedom play a cen- 
tral role. We show that the quarter-filled ground-state 
(GS) is a new insulating crystalline quantum state involv- 
ing resonant spin singlet electron pairs. We then argue, 
on the basis of numerical calculations, that pairing can 
emerge under light doping. Such a superconductor is ex- 
pected to break lattice symmetries in a similar way to the 
supersolid of Ref. [H (although with no charge ordering). 



One of the most exciting applications of these ideas is to 
assemblies of cold atoms in optical lattices, where models 
of the form we propose can in principle be realized by the 
appropriate tuning of atomic interactions 

The model: Our starting point is a (fermionic) Hub- 
bard model extended with nearest neighbor {{ij) in the 
sum) repulsion V on the checkerboard lattice: 



n = -tY,{flf,a + h. 

(1) 

using standard notation. We are interested in the limit 
where the on-site repulsion U is very large, forbidding 
double occupancy. Furthermore, we consider the limit 
V ^ t, in which case, for quarter filling ((n^) = i), the 
ground state is an insulator. Expanding about this state 
we obtain, to leading order, the effective Hamiltonian: 



□ 



where 



(2) 

L2 = ^ and the summation runs over the "empty 
squares" of the checkerboard lattice, with the sites of 
a given square counted ikjl. Ha acts on two electrons 
forming a singlet on the diagonal kl, transferring that 
singlet to the (empty) perpendicular diagonal ij. 

We also introduce a diagonal term which counts the 
squares where TYn can act 



= ( - - 2S, • S, ) n,nj{l - ni){l - n^) . (3) 



For W = t2 the Hamiltonian Ha + Hw becomes a sum of 
projectors, and we recover the physics of the Quantum 
Dimer Model at the Rokhsar-Kivelson point : the 
GS can be written exactly as an equal-weight superposi- 
tion of all zero-energy configurations [ld| . This nontrivial 
property also holds in the presence of static holes (t = 0). 
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FIG. 1: (Color online) Energies of the 1/4-filled GS and sin- 
glet excited states characterized by different quantum num- 
bers (as shown) vs W and J = (a). Inset: energy differ- 
ences w.r.t. GS. (b) Same data for excited states only and 
finite values of J. Phase transitions (shown by vertical bars) 
are signaled by the crossing between excited states (left) and 
by a W-independent energy (right). The dashed line in (a) 
shows the E^^^ perturbation result (see text). 



Finally, we also include the two-site spin exchange 



Hj 



Si • s 



(4) 



for [/ > y > t 



the exchange constant is given 
W + rrm-v^ ■ For the undoped 



where 

J — JJZTv 'W Uv U{U-V) 

system J and t2 are therefore independent parameters. 
Here we consider only antiferromagnetic (AF) coupling, 
J > 0, although a ferromagnetic coupling could also be 
realized for t <Q. 

Away from quarter filling, when Nh holes (or elec- 
trons) are introduced, simple counting shows that exactly 
2A'';i tetrahedra (named hereafter "half-hole tetrahedra" ) 
should contain a single, shared electron (or 3 shared elec- 
trons for electron doping). All other tetrahedra contain 
exactly two shared electrons, as before; violation of this 
constraint would lead to an energy increase of order V . 
The single hopping term t becomes then effective by mov- 
ing around the locations of these "fractional charges" 0] . 
We have performed extensive Lanczos Exact Diagonal- 
ization (ED) on a periodic (45° tilted) square cluster of 
A'^ = 32 sites in the insulator and for a small number of 
doped holes or electrons [ll| . 

Phase diagram at quarter-filling: For J = 0, by anal- 
ogy with Ref. [5|] and from simple analytic arguments, 
we expect to have as a function of W: (i) for large neg- 
ative W , parallel spin chains; (ii) for intermediate \W\ 
values a Resonating Singlet-Pair Crystal (RSPC); (iii) for 
W = t2,a. hquid-like "RK" point; (iv) for W > t2 a, man- 
ifold of isolated states (which include all ferromagnetic 
states), all having E = 0. This picture is indeed sup- 



ported by our ED calculations (see Fig. [TJa)) where the 
characterization of the various phases (see Fig. [2]) can be 
obtained from the analysis of the low-energy spectrum. 
In the two fold degenerate RSPC, electron pairs resonate 
in every second void plaquette, breaking translational 
symmetry. This is seen numerically in the collapse of a 
k = (tt, tt) Ai symmetry singlet excited state onto the 
GS [Fig.[Ha)]. For sufficiently negative W the electrons 
order along diagonal chains to optimize the Heisenberg 
exchange along the diagonals of the empty squares, hence 
breaking rotation symmetry. This is seen in the ED spec- 
trum as the quasi-degeneracy of the (k — 0,Bi) singlet 
and the GS. For W > 1, isolated states where the pla- 
quette resonance is ineffective (including the horizontal 
chains shown), are favored. 

Note that a perturbative expansion about isolated 
void plaquettes gives a very accurate estimate of the 
RSPC GS energy (for J — 0). Already in first order 
— £;(!) = N{^W - its), as shown in Fig. [T] — the 
agreement is excellent. A comparison to a similar expan- 
sion from disconnected diagonal chains gives an estimate 
of the boundary between the two phases in full agree- 
ment with the numerics, although it suggests that the 
actual critical value shifts to more neagtive W in the 
thermodynamic limit. One should point out that (i) in 
contrast to the bosonic case Q, the range of stability of 
the plaquette phase is much broader on the W axis and 
(ii) as shown in Fig. [TJb) and Fig. ^ the RSPC is also 
very stable when the exchange J is included. Only above 
J/t2 ^ 1.5, a new phase of 4-electron plaquettes oriented 
along the (1, ±1) directions is stabilized by a gain of ex- 
change energy (small Heisenberg loops have higher quan- 
tum fiuctuations) [l^. Lastly, we stress that our RSPC 
at J/t2 < 1.5 differs from the Bond Order Wave realized 
for t > at intermediate V/t ratio (for which charge fluc- 
tuations still occur) [l3| which shows two kinds of criss- 
crossed plaquettes. Indeed, although both phases exhibit 
a doubling of the unit cell (neither with charge ordering), 
they differ in the type of resonating plaquettes. 

Finite doping: In investigating the effect of doping, 
we restrict to = and consider only values of J for 
which the RSPC is realized at zero doping. 

A single doped hole (or electron), as discussed above, 
is in fact split into two mobile "half-hole tetrahedra" car- 
rying one electron only (and effective charge e/2) d, 0]- 
For reasons similar to those given in Ref. [6| , one expects 
a confining potential so that a small quasiparticle weight 
survives (in fact, in the t = Q limit, the two "half-hole 
tetrahedra" even share a common site). In order to esti- 
mate the coherent bandwidth we have computed by ED 
the GS energy of a single hole for various inequivalent 
k momenta |15| and results are reported in Fig. [3] For 
large t (compared to t2 or J) the bandwidth is reduced by, 
roughly, a factor t2/\t\ (or J/\t\) for t > 0, or even more 
for t < 0, and the hole becomes quite massive. This is 
very reminiscent of the behavior of a single hole doped in 
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FIG. 2: (Color online) Tentative phase diagram vs W and J. 
(See description of the singlet phases in the text). 
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FIG. 3: (Color online) Single hole bandwidth (in units of |t|) 
as a function of t'2.l\t\ for a fixed (physical) ratio of J jti. For 
comparison, the same quantity for the quantum AF on the 
square lattice is indicated, as a function of J/t, by a dotted 
line (blue online). 



a quantum antiferromagnet (shown also in Fig.[3]for com- 
parison) where the hole leaves behind a path of flipped 
spins which can be healed over a short time scale cx 1/J 
by spin fluctuations [l^ . We believe a similar path along 
which the plaquette order is perturbed exists also here 
behind the moving hole. 

As the eff'ective Hamiltonian lowers the energy of sin- 
glet pairs, it is natural to study the emergence of pairing 
interaction. For i = we found that two holes sit next 
to each other or, more precisely, that all four "half-hole 
tetrahedra" are fully packed around a single void plaque- 
tte in order to minimize plaquette resonance and bond 
exchange energies. The two-hole GS for t = can be very 
well approximated by — |vI/q^ where removes 

two electrons along the diagonals of a void plaquette with 
s-wave orbital symmetry. For finite t it is therefore conve- 
nient to define the overlap squared Z^h = |(^2h 
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FIG. 4: (Color online) Two hole Z-factor Z^h versus t'il\t\ for 
J = and J jti = 0.5. The case of two doped-electrons is also 
shown (.Z2e) as hashed symbols (for J = 0.5). For comparison 
the same quantity is shown as circles (blue online) for two 
holes doped into the quantum AF on the square lattice vs J/i 
(data for a 26-site cluster [18j). 



to get an accurate estimate of the so-called 2-hole Z- 
factor |(*2/i|A||^'o)p. Note that is now defined 

as a Bloch state with the same momentum as the 2 hole 
GS \^2h) (typically k = (0, 0)) i.e. as a linear superposi- 
tion of all local (degenerate) hole pair states. Increasing 
Z^h should remain finite as long as the holes stay 
bound. Fig. [5] shows that Z^h is more rapidly suppressed 
than for two holes in an AF [3l (shown also in Fig. [4] 
for comparison), suggesting that the hole pair is more 
fragile here for large |i|/i2 and |t|/ J values. A very sim- 
ilar behavior is also found for two doped electrons (see 
Fig. [3]) although the related bound state has now d-wave 
symmetry [l7|. 

To get a more direct insight on the tendency of a hole 
pair to break up at large |t|/i2 we have computed the 
pair binding energy A2/1 = Eih + -Bo/i - '^'Exh shown in 
Fig. O where Erah is the GS energy of the system with 
A^^ = m holes. From its definition, negative IS.2h suggests 
the binding of the holes, as previously seen in the doped 
Heisenberg AF on the square lattice (in blue) — even 
for rather large hole kinetic energy. On the contrary, the 
data for the 1 /4-filled checkerboard lattice do not provide 
evidence of binding when tjt^ > 1- Note however that 
A2/1 is subject to stronger finite size effects than Z2/1, so 
a weak pairing might still survive for \t\/t2 > 1. 

Since phase separation might compete with supercon- 
ductivity, one needs now to consider the four-hole bind- 
ing energy, A4 = E^h + Eoh — 2i?2/i, from which an esti- 

^r-§- at small hole density 

n-h — ^ can De given using k ""'^ ~ A'^A4. Phase separa- 
tion signaled by a negative curvature of E{nfi) implies, 
in this case, the formation of four-hole droplets (quar- 
tets) i.e. A4 < 0. Our numerical calculation of A4 on 
32 sites shows that quartets are stable when t is small 



mate of the compressibility k = 
— — can be given using k~ 
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FIG. 5: (Color online) Hole pair binding energy (in unit of t) 
vs t2/\t\ for both signs of t. For comparison, the case of the 
doped quantum AF on the square lattice is shown (data for 
N = 26) [il. 



enough, the plaquette resonance and bond exchange pro- 
viding an effective "glue" to bind holes together. Such 
a behavior is similar to the case of the doped quan- 
tum AF [l9| . However, here a very small hopping t can 
easily suppress phase separation while affecting pairing 
only slightly. Indeed, for a typical value of J/t2 — 0.5, 
A4 ~ — 0.0635 i2 (meaning k < 0) for static holes (t = 0) 
while A4 ~ -0.0042 12 (A4 ~ -0.0011 for t = 0.1 12 
{t = -0.1*2) and A4 ~ 2.471*2 (A4 ~ 2.044*2) for 
t — 0.5*2 (* — —0.5*2). We therefore expect a paired 
state when */*2 > 0.1. 

Lastly, we would like to draw possible scenarios at fi- 
nite doping: (i) If plaquette ordering survives at suffi- 
ciently small doping (plausible since the RSPC is "pro- 
tected" by a finite gap), one expects a superconduct- 
ing phase which inherits broken translational symmetry 
from its parent Mott insulator (the RSPC). This can be 
viewed as a new type of supersolid with no charge mod- 
ulation (20| : (ii) Alternatively, hole pairs could arrange 
along domain walls between out-of-phasc RSPC regions. 
This scenario is however unlikely, providing no obvious 
gain of kinetic energy w.r.t the supersolid. 

To summarize, on the checkerboard lattice, at quarter- 
filling (one electron per two sites) and in the limit of large 
NN repulsion (enforcing the "ice rule" constraint on ev- 
ery tetrahedron), the Hubbard model exhibits a robust 
insulating Resonating Singlet-Pair Crystal with a uni- 
form average charge. We provide arguments for Cooper 
pair formation and argue in favor of a supersolid phase in 
an extended region of parameters, phase separation being 
confined only at very small hole kinetic energy. However, 
hole-Cooper-pairs seem to be more fragile for large hop- 
ping than in the doped antiferromagnet on the square 
lattice. 
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